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Sunspots and the plethora of other phenomena occuring in the course of the 
11-year cycle of solar activity are a consequence of the emergence of magnetic 
flux at the solar surface. The observed orientations of bipolar sunspot groups 
imply that they originate from toroidal (azimuthally orientated) magnetic flux 
in the convective envelope of the Sun. We show that the net toroidal magnetic 
flux generated by differential rotation within a hemisphere of the convection 
zone is determined by the emerged magnetic flux at the solar surface and thus 
can be calculated from the observed magnetic held distribution. The main 
source of the toroidal flux is the roughly dipolar surface magnetic field at the 
polar caps, which peaks around the minima of the activity cycle. 

The basic concept for the large-scale solar dynamo involves a cycle during which the poloidal 
field and the toroidal field are mutually generated by one another (i, 2). The winding of the 
poloidal field by differential rotation creates a toroidal field. A reversed poloidal field results 
from the formation of magnetic loops in the toroidal field, which become twisted by the Cori¬ 
olis force due to solar rotation. In turn, the reversed poloidal field then becomes the source of 
a reversed toroidal field. In this way, the 11-year cycle of solar activity is connected to a 22-yr 
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cycle of magnetic polarity. 

Hale et al. (3) discovered that the magnetic orientations of the eastward and westward parts 
of bipolar sunspot groups in one solar hemisphere are the same during an 11-year cycle and 
opposite in the other hemisphere. This implies that the sunspot groups originate from a toroidal 
field of fixed orientation during a cycle. Toroidal flux of the opposite polarity would lead to 
sunspot groups violating Hale’s law. Since only a small minority of the sunspot groups are 
actually observed to violate this rule (4), opposite-polarity toroidal field is largely irrelevant as 
a source of sunspot groups. In other words, it is the hemispheric net toroidal magnetic flux 
given by the azimuthal average of the toroidal field that is relevant for the formation of sunspot 
groups. 

Here we use a simple method based on Stokes’ theorem to show that the emerged surface 
fields determine the net toroidal flux generated by differential rotation in a solar hemisphere. 
The time evolution of the net toroidal flux in the convection zone can thus be calculated using 
only observed quantities (differential rotation and field distribution at the surface). We compare 
the resulting net toroidal flux with the observed large-scale unsigned surface flux and find that 
they vary in a similar manner. 

We consider spherical polar coordinates, (r, 6*, 0), and the azimuthally averaged induction 
equation of magnetohydrodynamics, 

_ = V X (U X B +(u X b)-r/V X B), (1) 

where B(r, 0) and U(r, 0) are the 0-averaged magnetic field and plasma velocity, respectively, 
and rj is the magnetic diffusivity. Angular brackets indicate the azimuthal average. The term 
(u X b) denotes the correlation of the fluctuating quantities with respect to the azimuthal aver¬ 
ages, which gives rise to the a-effect and to enhanced (‘turbulent’) magnetic diffusivity (5). 

We define the contour 5T. enclosing the area S in a meridional plane of the Sun as shown in 
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Fig.m The direction of the contour is chosen such that the vectorial surface element of S points 
into the direction of positive azimuthal field, B^. Applying Stokes’ theorem to the integral of 
the induction equation over S yields the time derivative of the net toroidal flux, in the 
northern hemisphere of the convection zone, 

^ = = / (UxB + (uxb)-r 7 VxB)-dl, (2) 

at at V2s / j<5s 

where dS" is the surface element of S and dl is the line element along An analoguous 
procedure provides the net toroidal flux in the southern hemisphere, 

Rotation is by far the dominant component of the azimuthally averaged velocity, so that we 
write U = 0 = (fl r sin 9)^, where f2(r, 9) is the angular velocity and cj) the unit vector in 

the azimuthal direction. The effect of Xss(u x b) ■ dl reduces to that of the turbulent magnetic 
diffusivity, rjt, since the contribution of the a-effect to the generation of the toroidal field can be 
neglected against that of differential rotation (i). With we thus obtain 

d<f)^ f 

= (U X B - r/t V X B) ■ dl. (3) 

dr JsT, 

Guided by empirical results from helioseismology {6, 7), we take to be independent of r 
in the equatorial plane throughout the convection zone (5), i.e., f2(r, 7r/2) = fleq- This allows 
us to work in a reference frame rotating with angular velocity fleq, for which = 0 in the 
equatorial plane. We can further assume that the magnetic field does not penetrate the low- 
diffusivity radiative zone below the convection zone. Together with = 0 along the rotational 
axis, these assumptions imply that only the surface segment (d) in Fig. [T] contributes to the line 
integral of U x B along the contour We obtain 

/ (UxB)-dl= r^\sBrRQd9= f\Q-Qe„)BrRld{cos9), (4) 

JST. Jo Jo 

where U^, fl, and B^ are to be taken at the solar surface, r = Rq. This shows that the net toroidal 
flux generated in the convection zone by the action of differential rotation is determined by the 
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poloidal field threading the solar surfaee. Any additional poloidal flux that is fully eontained 
within the eonveetion zone would lead to equal amounts of East-West and West-East orientated 
toroidal flux, whieh do not eontribute to the net toroidal flux required by Hale’s law. 

The diffusion term in Eq. dH) is most relevant along the rotational axis, where toroidal flux 
ean be destroyed, and also in the equatorial segment of where flux ean be transported 
between the hemispheres. These proeesses, whieh deerease the net magnetie flux, are diffieult 
to quantify. In order to estimate their effeet, we approximate them with an exponential deeay 
term with e-folding time r. Altogether we then obtain 



(5) 


and an analoguous equation for the net toroidal flux generation in the southern hemisphere of 
the eonveetion zone. 

To evaluate the induetive part of Eq. ([5]) we use the synoptie magnetograms from the Kitt 
Peak National Observatory from 1975 to the present (9). The azimuthally averaged radial eom- 
ponent of the magnetie field, Br, as a funetion of time and cos 9 is shown in Eig. [2j4, where 
cos 9 = 0 represents the equator and cos 9 = ±1 the poles. The profile of the surfaee differen¬ 
tial rotation taken from {10), 


O — Oeq = —2.3 cos^ 9 — 1.62 cos^ 9 [°/day], 


( 6 ) 


given in Eig. [2^ is eombined with to obtain the quantitity (O — Oeq)-Br (Fig- [2t)> "whieh 
determines the generation of toroidal flux by differential rotation in Eq. ([5]). It is dominated by 
the eontribution from the polar regions. 

We ealeulated the time integral of the induetive part of Eq. ([5]) as well as its eounterpart 
for the southern hemisphere to obtain the net toroidal flux in both hemispheres as a funetion of 
time. The integration begins with zero flux in Eebruary 1975, the starting time of the synoptie 
observations. This is near solar aetivity minimum, during whieh time we expeet the toroidal flux 
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to change sign. The result given in Fig. [3] shows that the modulus of the toroidal flux generated 
from the polar fields reaehes peak values of the order of 1 — 6 x 10^^ Mx per hemisphere during 
recent activity cycles. Note that the net flux generated for the new eyele first has to eancel 
the opposite-polarity flux from the old eyele, so that it reaehes its peak value around activity 
maximum of the new eyele. 

The exponential deeay term in Eq. Q mainly leads to a phase shift of the time evolution: the 
sign reversals and the maximum values of the toroidal flux oeeur earlier sinee the flux from the 
previous cycle is eontinuously redueed by the decay. At the same time, the newly generated flux 
is also subjeet to the decay, so that the amplitude of the toroidal flux is only weakly affeeted: 
even for the extreme ease of r = 4 years, the peak values for cyeles 22 and 23 are redueed by at 
most 20-30%. Plots analoguous to Figure [3] for various values of r are provided in the online 
supplementary material. 

Is the inferred amount of net toroidal flux suffieient to provide the magnetic flux actually 
emerging in the form of bigger bipolar magnetic regions and sunspot groups at the solar surface? 
In order to estimate their eombined fluxes, whieh provide a rough measure of the subsurface 
toroidal flux, we ealeulated the integrated unsigned radial surface flux, $u = / |i?r|dA (where 
dA represents the solar surface element) for each hemisphere on the basis of the Kitt Peak syn- 
optie magnetograms. Because the spatial resolution element of the synoptie magnetograms is on 
average about 12 x 12 Mm^ on the solar surface, <hu represents the large-scale flux originating 
from the emergenee of bigger bipolar magnetie regions. Fig. |3] shows that the time evolution of 
•Fu for both hemispheres is very similar to that of the inferred toroidal flux, both with respect to 
phase as well as to amplitude (11). A short deeay time, r, would lead to a signifieant phase shift, 
suggesting that the deeay term in Eq. ([5]) has no strong effeet on the net flux. Fig.|3]indieates that 
the toroidal flux generated from the polar fields is suffieient to explain the observed time profile 
and amount of flux emerging in the resolved bipolar regions and sunspot groups in the course 
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of the solar cycle. Since there is probably considerable randomness in the flux emergence pro¬ 
cess, leading to elements of toroidal flux emerging more than once at various longitudes or not 
emerging at all, we cannot expect more than an order-of-magnitude agreement. 

Our results demonstrate that the emerged magnetic flux and particularly the polar fields are 
by far the dominating source of the net toroidal flux in the convection zone, from which the 
sunspot groups of the subsequent cycle originate. The solenoidality of the magnetic field means 
that the flux associated with the polar fields threads through the Sun’s convection zone, where 
it is wound up by differential rotation to generate toroidal magnetic flux. We thus confirm 
the conjecture of H.W. Babcock (72) that the observed polar fields represent the poloidal field 
source for the subsurface toroidal field. 

As also first suggested by Babcock, the net axial dipole moment represented by the po¬ 
lar fields results from the North-South dipole moments contributed by the individual sunspot 
groups and bipolar magnetic regions at the surface as a result of their systematic tilt with re¬ 
spect to the East-West direction (3). The tilt probably originates in one way or another from 
rotation via the Coriolis force: either by providing helicity to convective flows bringing mag¬ 
netic flux to the surface or by twisting buoyantly rising flux loops. The concept of Babcock was 
further developed by Leighton {13), who introduced the notion of surface flux transport for the 
buildup of the polar fields in connection with the dynamo process. More recently, surface flux 
transport models successfully reproduced the observed evolution of the surface fields and, in 
particular, the polar fields on the basis of the observed records of sunspot groups as flux input 
(14-18). This implies that the tilt of the larger bipolar magnetic regions determines the polar 
fields. Small bipolar regions and small-scale correlations are irrelevant in this respect. Together 
with the results shown here, this establishes the key part of the surface fields in the solar dynamo 
process and thus corroborate the basic dynamo concept of Babcock (72) and Leighton (13). 

The key role played by the polar fields for the generation of toroidal flux explains the strong 
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empirical correlation between the strength of the polar field {19) and the Sun’s open flux {20) 
around activity minimum with the number of sunspots of the subsequent activity cycle, which 
can be taken as a proxy for the underlying toroidal flux. Although the correlation is not perfect, 
which can be ascribed to randomness associated with the flux emergence process, it provides 
the best available method to predict the strength of the next cycle {21, 22). Our results put this 
method on a firm physical basis. 
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Figure 1: Integration eontour for the applieation of Stokes’ theorem. The eontour (thiek solid 
line) on a meridional plane of the Sun is used to ealeulate the net toroidal flux in the northern 
hemisphere generated by the aetion of differential rotation on the poloidal field. The thin solid 
line represents the solar surfaee, the dashed line the bottom of the eonveetion zone. The rotation 
poles are at the top and bottom. The eontour eonsists of a radial segment in the equatorial plane 
(a), a eireular are slightly below the bottom of the eonveetion zone (b), a part along the axis of 
rotation (e), and the solar surfaee (d). 
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Figure 2: A: Azimuthally averaged radial surfaee field from observed synoptie magnetograms 
as a funetion of cos 6 and time. B: Solar surfaee differential rotation relative to the equator as 
a funetion of cos 9. C: Map of the quantity (fl — representing the souree term for the 

generation of net toroidal flux in Eq. dH), as a funetion of cos 9 and time. 
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Figure 3: Observed magnetie flux at the solar surfaee and ealeulated net toroidal flux. Hemi- 
spherie total unsigned surfaee flux (dashed lines) from synoptie magnetograms and the modulus 
of the net toroidal flux (solid lines) are given for the last three solar eyeles. Red lines refer to 
the northern and blue lines to the southern hemisphere. 
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